Estimativa do erro de discretização analítico na solução de equações diferenciais utilizando o Método de Volumes Finitos

As análises de erros são conduzidas antes de qualquer projeto a ser desenvolvido. A necessidade do conhecimento do comportamento do erro numérico em malhas estruturadas e não-estruturadas surge com o aumento do uso destas malhas nos métodos de discretização. Desta forma, o objetivo deste trabalho foi criar uma metodologia para analisar os erros de discretização gerados através do truncamento na Série de Taylor, aplicados às equações de Poisson e de Advecção-Difusão estacionárias uni e bidimensionais, utilizando-se o Método de Volumes Finitos em malhas do tipo Voronoi. A escolha dessas equações se dá devido a sua grande utilização em testes de novos modelos matemáticos e função de interpolação. Foram usados os esquemas Central Difference Scheme (CDS) e Upwind Difference Scheme(UDS) nos termos advectivos. Verificou-se a influência do tipo de condição de contorno e a posição do ponto gerador do volume na solução numérica. Os resultados analíticos foram confrontados com resultados experimentais para dois tipos de malhas de Voronoi, uma malha cartesiana e outra triangular comprovando a influência da forma do volume finito na solução numérica obtida. Foi percebido no estudo que a discretização usando o esquema CDS tem erros menores do que a discretização usando o esquema UDS conforme literatura. Também se percebe a diferença nos erros em volumes vizinhos nas malhas triangulares o que faz com que não se tenha uma uniformidade nos gráficos dos erros estudados. Percebeu-se que as malhas cartesianas com nó no centróide do volume tem menor erro de discretização do que malhas triangulares. Mas o uso deste tipo de malha depende da geometria do problema estudado

Análise de erros da equação de advecção unidimensional no Método de Volumes Finitos

Uma análise utilizando a série de Taylor é apresentada para se estimar a priori os erros envolvidos na solução numérica da equação de advecção unidimensional com termo fonte, através do Método dos Volumes Finitos em uma malha do tipo uniforme e uma malha não uniforme. Também faz-se um estudo a posteriori para verificar a magnitude do erro de discretização e corroborar os resultados obtidos através da análise a priori. Por meio da técnica de solução manufaturada tem-se uma solução analítica para o problema, a qual facilita a análise dos resultados numéricos encontrados, e estuda-se ainda a influência das funções de interpolação UDS e CDS e do parâmetro u na solução numérica.

Desenvolvimento de um método de volumes finitos de alta ordem para a simulação de escoamentos de fluidos viscoelásticos

O uso da mecânica de fluidos computacional no estudo de processos envolvendo o escoamento de fluidos poliméricos está cada vez mais presente nas indústrias de transformação de polímeros. Um código computacional voltado a esta função, para que possa ser aplicado com sucesso, deve levar a predições mais próximas possível da realidade (modelagem), de uma forma relativamente rápida e eficiente (simulação). Em relação à etapa de modelagem, o ponto chave é a seleção de uma equação constitutiva que represente bem as características reológicas do fluido, dentre as diversas opções existentes. Para a etapa de simulação, ou seja, a resolução numérica das equações do modelo, existem diversas metodologias encontradas na literatura, cada qual com suas vantagens e desvantagens. Neste tópico se enquadra o trabalho em questão, que propõe uma nova metodologia para a resolução das equações governantes do escoamento de fluidos viscoelásticos. Esta se baseia no método dos volumes finitos, usando o arranjo co-localizado para as variáveis do problema, e na utilização de aproximações de alta ordem para os fluxos médios lineares e não-lineares e para outros termos não lineares que surgem da discretização das equações constitutivas. Nesta metodologia, trabalha-se com os valores médios das variáveis nos volumes durante todo o processo de resolução, sendo que os valores pontuais são obtidos ao final do procedimento via deconvolução. A solução do sistema de equações não lineares, resultante da discretização das equações, é feita de forma simultânea, usando o método de Newton São mostrados então, resultados da aplicação da metodologia proposta em problemas envolvendo escoamentos de fluidos newtonianos e fluidos viscoelásticos. Para descrever o comportamento reológico destes últimos, são usadas duas equações constitutivas, que são o modelo de Oldroyd-B e o modelo de Phan-Thien-Tanner Simplificado. Por estes resultados pode-se ver que a metodologia é muito promissora, apresentando algumas vantagens frente às metodologias convencionais em volumes finitos. A implementação atual da metodologia desenvolvida está restrita a malhas uniformes e, consequentemente, soluções para problemas com geometrias complexas, que necessitam de refinamento localizado da malha, foram obtidas somente para baixos números de Weissenberg, devido a limitação do custo computacional. Esta restrição pode ser contornada, tornando o seu uso competitivo.

Simulação numérica de escoamentos de fluidos pelo método de elementos finitos baseado em volumes de controle em malhas não estruturadas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Simulação numérica de escoamentos de fluídos pelo método de elementos finitos baseado em volumes de controle com a técnica de passo fracionado

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Simulação de grandes escalas de escoamentos incompressíveis com transferência de calor e massa por um método de elementos finitos de subdomínio

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Métodos numéricos para ecuaciones diferenciales rígidas. Aplicación a la semidiscretización del método de elementos finitos

295 p.

Análise da dissipação das tensões em dentes humanos restaurados com facetas laminadas de cerâmica, com três tipos de preparos, através do método dos elementos finitos

O objetivo deste trabalho é analisar in vitro a dissipação de tensões em incisivos centrais superiores humanos restaurados com facetas de cerâmica feldspática, através da análise do método dos elementos finitos, considerando cargas funcionais de mastigação e corte dos alimentos, em função de três tipos de preparos utilizados: sem proteção incisal; com proteção incisal em ângulo e com proteção incisal em degrau palatino. Foram utilizadas modelagens bidimensionais de um incisivo central superior e suas estruturas de suporte, simulando três situações: (Primeira modelagem) incisivo central superior com desgaste vestibular (em forma de janela); (Segunda modelagem) incisivo central superior com desgaste vestibular e proteção incisal em plano inclinado; (Terceira modelagem) incisivo central superior com desgaste vestibular, e proteção incisal com degrau palatino. Foi considerada uma carga (P=100N) com uma inclinação de 45 concentrada, simulando a região de contato do incisivo central inferior com o superior durante a mastigação e uma na região de contato topo a topo dos incisivos superior e inferior, simulando o corte dos alimentos. Após a análise dos dados obtidos pela distribuição de tensões, pode-se concluir que quanto à dissipação das tensões em todo o sistema proposto, com a aplicação de carga em 45, não foram observadas mudanças no estado tensional nos três diferentes preparos. Quando foi aplicada carga vertical, simulando o contato de topo, houve variação no estado tensional no sistema do dente com preparo em janela. Nas facetas, com a aplicação de carga em 45, nos preparos em janela e com proteção incisal em plano inclinado o resultado foi semelhante nos valores tensionais enquanto, nas facetas em dentes preparados com proteção incisal com degrau palatino, a distribuição foi mais homogênea tendo valores superiores, mostrando que o abraçamento do dente diminuiu a flexão.

Avaliação de ligações soldadas de perfis tubulares em estruturas de aço através do método dos elementos finitos.

A necessidade de se obter soluções para os diversos problemas estruturais na engenharia, associado ao amplo crescimento e difusão da construção metálica fez surgir, por volta dos anos 60, os perfis tubulares. As ligações soldadas entre perfis tubulares são as ligações de maior difusão nas construções com este tipo de perfil, sendo necessário portanto, estudos mais detalhados que tornem a utilização de ligações soldadas mais seguras e otimizadas. O Eurocode 3 possui um item específico para as ligações entre perfis tubulares, em especial as ligações soldadas devido a suas particularidades. Infelizmente, a norma brasileira NBR8800 não aborda as ligações envolvendo perfis tubulares. Para a execução da análise numérica no presente trabalho, faz-se necessário uma revisão bibliográfica. Com estes resultados, foram feitas modelagens de ligações tubulares com o método dos elementos finitos de forma a otimizar os modelos a serem utilizados em uma análise paramétrica futura. Desta forma, foram desenvolvidos dois modelos numéricos, um considerando ligação tipo T entre perfis tubulares quadrados e outro para uma ligação tipo K entre perfis circulares. Estes modelos foram caracterizados através de elementos de casca com seis graus de liberdade por nó considerando-se adicionalmente o efeito de membrana. A análise não-linear realizada considerou a não-linearidade do material através do critério de plastificação de Von Mises através de uma lei constitutiva tensão versus deformação bi-linear e a não-linearidade geométrica foi atribuída através da formulação de Lagrange atualizado. Dentre as principais conclusões obtidas no presente trabalho, pode-se citar que os resultados para as ligações tipo T, o Eurocode 3, fornece resultados que precisam ser observados com cautela. Todavia, para as ligações do tipo K, os resultados numéricos mostraram-se sempre inferiores aos valores através do Eurocode 3, representando um dimensionamento a favor da segurança.

Desenvolvimento e implementação de algoritmo para análise de dutos fabricados ou reparados com material compósito carregados com pressão interna utilizando o método dos elementos finitos.

A integridade de dutos ganha importância à medida em que o desenvolvimento da indústria expande a malha para transportar líquidos e gases num contexto de aumento das exigências ambientais. Um importante aliado para manutenção da integridade de dutos são reparos de materiais compósitos. Estes materiais apresentam baixa densidade, capacidade de direcionar resistência de acordo com as direções das solicitações, execução de reparo a frio sem necessidade de interromper produção ou grande maquinário. Este trabalho tem como objetivo desenvolver e implementar um algoritmo de elementos finitos que permita avaliar os esforços e a resistência das paredes de um tubos fabricados ou reparados com laminados de material compósito carregados com pressão interna. Entre as vantagens de desenvolver um programa tem-se: agilidade de avaliação, menor custo com licença, menores exigências computacionais, possibilidade de desenvolver o programa e o melhor entendimento da modelagem dos fenômenos. Utiliza-se como entrada do programa o diâmetro do duto, pressão interna e parâmetros do laminado. A modelagem em elementos finitos é realizada a partir da teoria clássica de laminados. Aplicando o carregamento resultante da pressão interna, determina-se os deslocamentos e são calculadas as tensões e aplicado o critério de falha de Tsai-Hill em cada camada. Estudos experimentais e numéricos encontrados na literatura foram simulados com o programa gerado e os resultados para propriedades do laminado, tensões nos dutos e pressão de ruptura apresentam concordância com os resultados da literatura.O programa ainda tem sua estrutura modificada para encontrar a pressão de falha a partir dos dados do laminado. O programa implementado permite uma avaliação rápida de resistência do reparo e possibilita avaliar rapidamente a resposta a mudanças nos parâmetros de projeto do laminado.

Projecto de equipamento de rotação de moldes de grande porte com recurso ao método de elementos finitos

A fabricação de moldes é uma actividade que se desenvolve em Portugal há várias décadas, com qualidade reconhecida mundialmente. A indústria automóvel é um dos principais destinos de grande parte dos moldes produzidos em Portugal. Algumas das peças produzidas nesses moldes possuem dimensões apreciáveis, tendo em conta o peso e dimensão média das peças injectadas em plástico. São exemplo disso os pára-choques dos automóveis e o tabliê. A produção dos moldes para estas peças é complexa, implicando um longo tempo de maquinagem e a manipulação dos moldes em diferentes fases, com vista ao acesso a todas as faces do molde. A IGM – Indústria Global de Moldes, SA. é a empresa responsável pela produção dos moldes para peças de média dimensão dentro do Grupo SIMOLDES. Atendendo à necessidade permanente de rodar a 90o moldes de grande porte, que podem apresentar pesos na ordem das 10 a 30 toneladas, e não existindo no mercado qualquer solução adaptável a esta necessidade, a empresa entendeu por bem levar a efeito o projecto desse equipamento, atendendo ao compromisso custo-benefício que torne viável a realização prática do mesmo. Após os esboços iniciais e uma discussão interactiva e iterativa com a empresa, foram analisadas as diferentes soluções entendidas como viáveis, tendo sido escolhido um dos anteprojectos realizados. Foram ainda discutidas as diversas alternativas de accionamento. Com base nesse anteprojecto, a estrutura foi optimizada e verificada através do Método de Elementos Finitos, tendo sido elaborado o projecto final, com o grau de detalhe necessário à sua fabricação e implementação na empresa.

Criação de programa informático em ambiente Matlab® para estudo de estruturas planas pelo Método de Elementos Finitos.

Os modelos a ser analisados pelo Método de Elementos Finitos são cada vez mais complexos e, nos tempos que correm, seria impensável realizar tais análises sem um apoio computorizado. Existe para esta finalidade uma vasta gama de programas que permitem realizar tarefas que passam pelo desenho de estruturas, análise estática de cargas, análise dinâmica e vibrações, visualização do comportamento físico (deformações) em tempo real, que permitem a otimização da estrutura. Sob o pretexto de permitir a qualquer utilizador uma análise de estruturas simples com o Método dos Elementos Finitos, surge esta tese, onde se irá criar de raiz um programa com interface gráfica no ambiente MATLAB® para análise de estruturas simples com dois tipos de elemento finito, triangular de deformação constante e quadrangular de deformação linear. O software desenvolvido, verificado por comparação com um software comercial dedicado para o efeito, efetua malhagem com elementos bidimensionais triangulares e quadriláteros e resolve modelos arbitrados pelo Método de Elementos Finitos, representando estes resultados visualmente e em formato tabular.

Determinação das tensões residuais e deformações resultantes do processo de soldadura TIG através do método dos elementos finitos

Dissertação para obtenção do Grau de Mestre em Engenharia Mecânica

Aplicaciones del método de elementos finitos a problemas de termofluídos

Tesis (Maestría en Ciencias de la Ingeniería Mecánica con Especialidad en Térmica y Fluídos) UANL

Simulación del proceso de forja de codos sin costura y evolución de la herramienta de formado por el método de elementos finitos

Tesis ( Maestro en Ciencias de la Ingeniería Mecánica con Especialidad en Materiales ) U.A.N.L.